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Abstract 

We present a one-loop calculation of the static potential in the SU(2)-Higgs model. The 
| connection to the coupling constant definition used in lattice simulations is clarified. The 

I ' consequences in comparing lattice simulations and perturbative results for finite temperature 

Q_i' applications are explored. 

m : 

^ '. 1 Introduction 

i ' One of the most fundamental and intriguing problems in present-day particle physics is the observed 
^ ■ baryon asymmetry of the universe — that is the fact that we do not see hardly any antimatter around 
I/* - ) . us, but we can see a large amount of matter, and a vast amount of photons. 

One could be inclined to attribute this asymmetry to the initial conditions of the big bang — 
q [ but this being rather ad hoc, we shall focus on dynamical generation of baryons. Another possibility 
0\ could be arguing that the entire universe has no net baryon number, just our region is more baryonic, 
^ ■ while other regions are more antibaryonic. However, from the boundary between two such regions, 
^ ■ characteristic photons would bring us information, and lacking this, we can safely state that in our 
■ ! vicinity of at least 10 13 solar masses there is just matter. Since no known mechanism can separate 
<D \ matter and antimatter on such large scales, we shall assume that the universe is baryonic and this 
^ ' asymmetry was formed some time after the big bang. 

In 1967 Sakharov gave three necessary conditions for baryogenesis, which are 



X 



1. Baryon number violation — which is obvious 

2. C and CP violation — otherwise the number of baryons and antibaryons generated would be 
equal 

3. Departure from thermal equilibrium — since quantum numbers do not change in thermal equi- 
librium 

Baryon number violating processes are known to be present in GUT's — and GUT's are still favoured 
candidates for baryogenesis. However, it came as a slight surprise that even within the standard 
model there are such processes, the nonperturbative sphaleron processes. The standard model be- 
ing experimentally verified with great precision, the investigation of the possibility of electroweak 
baryogenesis is clearly of great importance. 

The departure from thermal equilibrium can be realised during the electroweak phase transition. 
This is customarily described by the effective potential. 
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Figure V. The effective potential as a function of the scalar field. The phase transition point is 
defined by the degeneracy of the two minima 



In order to have electroweak baryogenesis, the phase transition must be strongly first order, or 
more quantitatively the relation 

m/T C > i 

must hold, i.e. the vev of the Higgs field has to exceed the critical temperature. 

How can we study this relation? The most straightforward method is resummed perturbation 
theory (cf. e.g. [I], ||, U). In the low temperature Higgs phase the perturbative approach is expected to 
work well, however, serious infrared problems are present in the high temperature symmetric phase. 
Since the determination of thermodynamical quantities at the critical temperatures is based on the 
properties of both phases, non-perturbative techniques are necessary for a quantitative understanding 
of the phase transition. 

One very successful possibility is to construct an effective 3-dimensional theory by using di- 
mensional reduction, which is a perturbative step. The non-perturbative study is carried out in this 
effective 3-dimensional model (see e.g. |4| and references therein). Analytical estimates are confirmed 
by numerical results and relative errors are believed to be at the percent level. 

Another approach is to use 4-dimensional simulations. The complete lattice analysis of the 
standard model is not feasible due to the presence of chiral fermions. however, the infrared problems 
are connected only with the bosonic sector. These are the reasons why the problem is usually studied 
by simulating the SU(2)-Higgs model on 4-dimensional lattices, and perturbative steps are used to 
include the U(l) gauge group and the fermions. 

Despite the fact that both perturbative and lattice approaches are systematic and well-defined, it 
is not easy to compare their predictions. The reason for this is that in lattice simulations the gauge 
coupling constant is determined from the static potential, whereas in perturbation theory the MS 
scheme is used. Therefore, if one wishes to compare perturbative and lattice results, a perturbative 
calculation of the static potential proves to be of great help. 
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I Calculation or trie one-loop static potential 

The concept of the static potential was introduced more than 20 years ago by L. Susskind 0. 



Figure 2: Heavy quark-antiquark loop 



When a very heavy quark-antiquark pair is created from vacuum, separated and kept at distance 
r for time t and then let annihilate, the matrix element of the process can be given as 

(i\e- m \f) = e- v ^(i\f), 

and so 

7(r) = - + lim 7 log(W P , t ). (1) 

The Feynman rules for these static sources are quite simple. Owing to their great mass, their 
coordinate space propagator is purely timelike, 



pT 

b 

i S Q (p) i S A (p) 



iS Q / A {x -y) = 6(x- y)0(±zo =F Vo), 



which upon Fourier transformation gives the momentum space propagator 



iS$ /A (p) 



vp + ie 



where v is the velocity of the source ([1,0,0,0] to a first approximation). 
The heavy ( ant i) quark-gauge boson vertex is given by 
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By means of these rules, the one-loop static potential was calculated long ago in quantum elec- 
trodynamics and quantum chromodynamics [||, |6], [7| , and even the full two-loop result was published 
recently ||. The case of QED is rather simple: it was shown in [|] that summing up all orders in 
perturbation theory is equivalent to taking the exponential of the one gauge boson exchange graph. 
However, this calculation is more than just a warm-up excercise, as the abelian parts of the graphs 
in more complicated theories (QCD, SU(2)-Higgs model) are identical to those of QED. Therefore 
we can focus on graphs with nonabelian contributions. 

Our calculation was performed in the MS scheme and the Feynman gauge but the result is gauge 
independent, as it should be for a physical observable. The relevant graphs are shown in Fig. 3. 

Solid lines represent the heavy quark (antiquark) propagator, while wavy lines the vector boson 
propagator. External heavy quark (antiquark) propagators are not shown in the figure. The one-loop 
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Figure 3: Graphs giving nonvanishing contributions to the static potential 
corrected vector boson propagator contains scalar and ghost contributions as well. The result can 



be conveniently given in momentum space. One obtains || [TO 
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where k 2 denotes the square of the three-momentum k, Mu the Higgs mass and Rhw 
The function F is defined as 
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Although the dependence on the renormaization scale can be removed by introducing the one-loop 
W mass |10|| , we do not follow this line. 

Eq. (Q) has to be Fourier transformed into coordinate space. We applied brute force methods 
performing numerical integration. As a check, we compared our results with various pieces of the 



partly analytic calculation in |T0[ for the derivative of the potential (with respect to distance). The 
agreement is excellent. 

Our result is presented in Figs. 2 and 3, where the various parts of the one-loop correction to the 
potential are plotted. We define 



4 



B 



r*M» 



Figure 4: The coefficients of g 4 /(167r 2 ) — curve A 
static potencial Eq. (B) as a function of distance times W mass. 
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where M$/ = My/ — 5 My/, with 5Mvi/ the one-loop mass correction. Since 5 My/ is scale dependent, 
so is Myj. A and B are functions of the distance r and Rhw = Mh/Mw- We choose My/ = 80GeV. 
Fig. 4 shows the dependence of A and B on the dimensionless distance rMy/ for Rhw = 0.8314 
(corresponding to the end point of the first order finite temperature phase transition JlTJ), while 
Fig. 5 shows the Rhw dependence for r = M w . 



3 Relation of the continuum version of the lattice coupling 
constant definition to the MS coupling constant 

Since we wish to compare results of lattice simulations and continuum perturbation theory calcu- 
lations, it is an essential point to define the SU(2) gauge coupling in the same way in both cases. 
However, in continuum perturbation theory the MS running coupling constant at a given renormal- 
ization scale is more natural (as used in Eqs. (||]}4|), too), while in lattice simulations other definitions 
are applied. Therefore we have to establish the relation between the coupling constants. 

The lattice definition of the coupling constant (inspired by [|T^]) is given in fl~3|| . 

First rectangular Wilson loops of size (r, t) are measured. Extrapolating to large t and dividing 
the logarithm by —t one gets the static potential in the t — > oo limit as a function of r [see Eq. ([[])]. 
The nonperturbative lattice static potential is fitted by a finite lattice version of the Yukawa potential 
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Figure 5: The coefficients of g 4 /(167r 2 ) — curve A — and of g 4 /(16n 2 ) log(/i 2 /M^/) — curve B — of the 
static potencial Eq. (|) as a function of Rrw- The distance is M^} 1 . 



with four parameters (for details cf. One of these parameters is the mass in the exponential of 

the Yukawa potential, which is usually called the screening mass. The gauge coupling at distance r 
is defined as the ratio of the discrete r derivative of the lattice simulated nonperturbative potential 
and the discrete derivative of the tree-level lattice Yukawa potential normalized by the square of 
the tree-level coupling and with the mass parameter Mi attice identified with the screening mass. In 
practice g 2 attice (lattice) i s determined and is called the local renormalized gauge coupling constant 
on the lattice. The lattice results at various Higgs masses are collected in Table 1. Data are from 

ni, m 0, m ■ 



The gauge coupling constant can be defined in the same spirit in the case of continuum pertur- 
bation theory, too: 

i 4- {- v ^ 

9 2 r(t) = —~ ^ T^T, (5) 
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i.e. by taking the ratio of the derivatives with respect to r of the one-loop potential and the tree-level 
potential normalized by the square of the tree-level coupling. In Eq. ( 
Cp = 3/4, and M screen is obtained from the fit. 

Since M screen — = 0(g 2 ), for distances satisfying M screen — 1/r 
into the form 



V(r) 
0(g 2 ) 



is given by Eq. 
we can put Eq. 
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C and D are functions of Rhw and M screen , 
80GeV. 



their values are tabulated in Table 2 for M s , 
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Table 1: Various quantities calculated for values of Rhw used in lattice simulations. 



In this procedure we have to choose the gauge coupling in the one-loop potential so that 
^(Mg^een) reproduces the lattice result (third row of Table 1) for the appropriate value of the 
Higgs mass. For the applications of the following section (thermodynamical quantities at and around 
the critical temperature Tc of the first order electroweak phase transition) the scale of the one-loop 
potential is chosen to be Tc ~ 2M#, where Mh is the Higgs boson mass at zero temperature. Thus 
the gauge coupling appearing in the one-loop potential is actually the MS gauge coupling at scale 
Tc- The MS gauge coupling values obtained from this procedure are given in the sixth row of Table 
1. 

Other definitions of the perturbative gauge coupling are also possible [10]; however, our definition 
seems to provide the most systematic way of comparing perturbative and nonperturbative results. 
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Table 2: The coefficients C and D defined in Eq. (|j) as a function of Rhw- 



4 Comparison of perturbative and lattice results for physi- 
cal observables 

In this section we compare lattice results and perturbative predictions for the finite temperature 
electroweak phase transition. 
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Lattice Monte Carlo simulations provide a well-defined and systematic approach to study the 
features of the finite temperature electroweak phase transition. During the last years large scale 
numerical simulations have been carried out in four dimensions in order to clarify non-perturbative 



details |I3|,[n}|,||I5|,|fR)|. Thermodynamical quantities (e.g. critical temperature, jump of the order 
parameter, interface tension, latent heat) have been determined and extrapolation to the continuum 
limit has been performed in several cases. Nevertheless, it has proven difficult to compare perturbative 
and lattice results, because the perturbative approach used the MS scheme for the gauge coupling, 
whereas the lattice determination of the gauge coupling has been based on the static potential. This 
difference between the definitions can be removed on the basis of the previous section. 

In this paper we use the published perturbative two-loop result for the finite temperature ef- 
fective potential of the SU(2)-Higgs model ||. Note that the numerical evaluation of the one-loop 
temperature integrals gives a result which agrees with the approximation based on high temperature 
expansion within a few percent. The reason for this is that the perturbative expansion up to order 
g A , A 2 corresponds to a high temperature expansion, which is quite precise for the Higgs boson masses 
we studied. It is known that the perturbative loop expansion becomes unreliable for Higgs masses 
above approximately 50 GeV (e.g. resummed perturbation theory fails to predict the end-point of 
the electroweak phase transition, thus it gives a first order phase transition for arbitrarily large Higgs 
boson masses). In the physically relevant range of the parameter space the electroweak phase tran- 
sition can only be understood by means of non-perturbative methods. Therefore it is particularly 
instructive to see quantitatively how perturbative and lattice results agree for small Higgs boson 
masses and how they differ for larger ones. 

In lattice simulations masses are extracted from correlation functions, and it is possible to use 
the zero temperature effective potential in order to include the most important mass renormalization 
effects. The Higgs boson mass obtained from the asymptotics of the correlation function corresponds 
to the physical mass determined by the pole of the propagators, i.e. the solution of p 2 — M 2 = U(p 2 ), 
where n(p 2 ) is the self-energy. The effective potential approach suggested by Arnold and Espinosa 
approximates n(p 2 ) by n(0) in the above dispersion relation. It has been argued that the difference 
between the two expressions is of order g 5 v 2 (v is the zero-temperature vacuum expectation value), 
which does not affect our discussion. In this scheme the correction to the MS potential reads 

'dm 2 + -L$a) + %\ (7) 





where 

256vr 2 256vr 2 \ fi 

Here /i is the renormalization scale and My/ is the W-boson mass at T = 0. The above notation 
corresponds to a tree-level potential of the form m 2 ip 2 /2 + \(p A /4. Note that this treatment is 
analogous to previous comparisons of perturbative and lattice results Jl6| . 

In [11, 0, several observables were determined, including renormalized masses at zero 



temperature (M H , M w ), critical temperatures (T c ), jumps of the order parameter (y?+), latent heats 
(Q) and surface tensions (a) for different Higgs boson masses. As usual, the dimensionful quantities 
were normalized by the proper power of the critical temperature. Simulations were performed on 
L t = 2, 3, 4, 5 lattices (L t is the temporal extension of the finite-temperature lattice) and whenever 
it was possible a systematic continuum limit extrapolation was carried out assuming standard 1/a 2 
corrections for the bosonic theory. 

The statistical errors of these observables are normally determined by comparing statistically 
independent samples. Jackknife and bootstrap techniques were used |17| and correlated fits were 



performed |T3| to obtain reliable estimates of the statistical uncertainties. A correct comparison has 
to include errors on the parameters used in the perturbative calculation. These uncertainties are 
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Table 3: Comparison of the perturbative and the lattice results. 



connected with the fact that neither the Higgs boson mass nor the gauge coupling constant can be 
determined exactly in lattice simulations. Including these errors, the perturbative prediction for an 
observable is rather an interval than one definite value. 

To obtain a better measure of the correspondence between perturbative and nonperturbative 
results, and to incorporate their errors, one introduces "pulls" defined by the expression 

perturbative mean — nonperturbative mean . , 

Pull = : : : : • 9 

perturbative error + nonperturbative error 

The four different pulls at different Higgs boson masses are tabulated in Table 4 and plotted in Fig. 6. 
For the sake of convenience, we used the shorthand P T = pull of T C /M H , P4, = pull of y? + /T c , 
P Q = pull of Q/T£, and P a = pull of a/T%. 



m H (GeV) 


16.4(7) 


33.7(10) 
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Table 4: Values of the four different pulls for various Higgs boson masses 



The quantity which has the smallest pull even for large Higgs boson masses is T C /M H - A quadratic 
fit was performed to this quantity as a function of Rhw- The result is 



~ C - 2.494 - 0M2R HW + 0.223R 2 HW . (10) 



M: 



H 



For large Higgs masses he unreliability of perturbative predictions (in particular concerning the 
quantity P^, which is of central importance from the viewpoint of baryogenesis) is striking. 



5 Conclusions 

Searching for the origin of the observed baryon asymmetry of the universe one should focus on the 
electroweak phase transition. By calculating the SU(2)-Higgs static potential perturbatively, one can 
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Figure 6: "Pulls" plotted against the Higgs mass. Arrows indicate values outside the interval [—5, 5]. 



establish a better connection between perturbative and nonperturbative studies of the electroweak 
phase transition. 

From this relation it can be seen that a purely perturbative study is not satisfactory. One also 
comes to the conclusion that electroweak baryogenesis is ruled out in the standard model — or to put 
it rather positively: there is physics beyond the standard model. 

Therefore, extensions of the standard model have to be studied — preferrably on the lattice, either 
in dimensionally reduced theories or in 4D. A very promising candidate is the MSSM, although the 
large number of free parameters makes this study a formidable task. However, within the MSSM 
the baryon asymmetry can be accounted for only if the lightest Higgs mass is Mh < 115GeV (for a 
|19|). Therefore experiments will either rule out this scenario in the close future, 



recent review, see 



or, measuring several parameters, will facilitate numerical simulations. 
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